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1. Scaling Laws for Fluxon Production 

1.1. BACKGROUND 

Although equilibrium, or adiabatic, correlation lengths £ a d(T) diverge at 
the critical temperature T c of continuous phase transitions, correlation 
lengths always remain bounded, in practice. This is because causality pre- 
vents a system becoming ordered on very large scales within the finite time 
in which transitions are implemented. In consequence, the order parameter 
fields become frustrated, and defects arise to mediate between the different 
equivalent ground states of the system. By observing these defects (in our 
case, fluxons) we obtain a direct experimental guide to the way in which 
the transition has been implemented. 

Causal bounds on defect densities produced in the early universe at 
GUT transitions were proposed by Kibble[l, 2] but, because of uncertainty 
in the models, these predictions are not robust. However, Zurek suggested[3, 
4] that identical causal arguments were applicable to condensed matter 
systems, for which direct experiments on defects could be performed. 

The reasoning is that, as the temperature T(t) falls through T c (at time 
t = 0, say), there is a maximum speed c(t) = c(T(t)) at which the sys- 
tem can become ordered. Causality suggests that a domain structure, with 
defects linking domains, cannot be established before time i, determined 
by |£ad(£)| = c(t), the first time at which the rate of collapse of the adia- 
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batic correlation length is comparable to the maximum permitted speed. 
For times t < i the system is frozen, and t measures the first time that we 
re-enter an adiabatic regimefl, 2, 3, 4]. In appropriate units t = ^{tq/tqY , 
where the exponent v depends on the system, and tq is the quench time. 

The second ingredient of the Kibble-Zurek analysis was to identify the 
defect separation at their time of formation as 



The value of a depends on the system, as do £o and To- It is this scaling 
behaviour that we have tested in our experiment, using annular Josephson 
tunnel junctions (JTJs), for which the defects are fluxons. 

In fact, the situation is more complex. Defects can only appear with 
the correct energy profiles when the order parameter has achieved its final 
state magnitude. This is effected by the exponential growth of unstable 
modes and is largely completed by the spinodal time t sp , for a transition 
implemented by phase separation. However, for simple systems it can be 
shown [5, 6] that, although t sp > t, t sp = t, up to 0(1) multiplicative 
logarithmic corrections in tq and the microscopic parameters of the system. 

Further, for simple cases, when defects are sheets, tubes or balls of the 
unstable ground state, the separation of defects is not governed directly 
by £,ad, derived from the large-distance behaviour of the field correlator. 
Rather, it is governed by the separation £ 2ero of the field zeroes that label 
these false groundstates, determined by the short-distance behaviour of the 
correlation functions. We shall see that this very different understanding of 
how defects form is crucial for our analysis, when we consider systems of 
linear size L < £. Nonetheless, the differences between (, a d(t) and (, Z ero(t S p) 
still remain logarithmic in tq, and the simple scaling behaviour (1) is reli- 
able at the level of accuracy of the experiments [7, 8, 9, 10, 11, 12, 13] that, 
prior to ours, have been performed to check it. 

1.2. THE SCALING PREDICTIONS FOR FLUXONS 

The order parameter of a Josephson tunnel junction at temperature T < 
T c is the phase difference <p of the macroscopic superconducting quantum 
mechanical wave functions across the barrier. For an annular JTJ (AJTJ) 
with a bias current T, (ft obeys [14] the perturbed Sine-Gordon equation 
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provided the width Ar of the annulus, of radius r, satisfies Ar r and 
Ar Aj(T), the Josephson coherence length; x measures the distance 
along the annulus, c(T) is the Swihart velocity, and a and j3 are the 
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coefficients of the losses due to the tunneling current and to the surface 
impedance, respectively. 

The boundary conditions for Eq.(2) are periodic, <p(x + C) = (f>{x) + 2im, 
where C = 2irr is the circumference of the junction and the winding number 
n is an integer corresponding to the algebraic sum of fluxons trapped in 
the junction barrier at the normal-superconducting (N-S) transition. 

The classical fluxons are the 'kinks' of the Sine-Gordon theory. Eq. (2) 
is only valid once the transition is complete, and we shall not use it to study 
the appearance of fluxons. However, it is sufficient to enable us to identify 
Aj(T), diverging at T c , as the equilibrium correlation length £ a d(T) to be 
constrained by causality. Further, the Swihart velocity c(T) measures the 
maximum speed at which the order parameter can change[15, 16]. 

For the experiment that we shall describe below, we find that C < £. In 
fact, because we are really counting field zeroes (mod 2ir), the conclusions 
are equally valid for small systems. Any field crossing zero (mod 2tt) has 
the potential to mature into a fluxon. Thus, before the causal time we 
have a picture in which there is a fractal thermal fuzz of potential fluxons, 
whose density depends on the scale at which we look. By 0(F) some of 
these have developed into the fluxons that we see subsequently (see [5, 6]). 
However, because from this viewpoint we have many proto-defects jockeying 
to become the real thing, the relevant scale to compare to the system size 
is not £ but £o> and C S> £o- 

A detailed discussion of JTJs has been given elsewhere by us [17, 18], 
including a summary of this experiment [19], and we refer the reader to 
these articles for more details. The JTJs in our experiment are symmetric, 

1. e. the electrodes are made of identical superconducting material. For such 
JTJs v = 1/2 and a = 1/4. Therefore, at the time of their formation the 
separation of fluxons is expected to be 

^(tq/to) 1 / 4 . (3) 

In terms of the parameters of the JTJs, £o = \Zh/2efiod s aJ c (0), where 
J C (T) is the critical Josephson current at temperature T for a JTJ with 
superconductors with effective thickness d s . The parameter a is given in 
terms of the gap energy and critical temperature and has a value between 
3 and 5. As for tq, it is given as To = £o/co> where Co defines the behavior 
c(t) = co(t/rQ) 1 / 2 of the Swihart velocity for the system near T = T c . 

2. The Experiment 

2.1. MEASURING FLUXONS 

When fluxons are produced they are static. To see them it is necessary to 
feed a bias current to the AJTJ, whereupon they move as magnetic dipoles. 
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Figure 1. The measured current- voltage characteristics of an AJTJ without trapped 
fluxons. For each current branch the number of fluxon-antifluxon pairs FF is indicated. 
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Figure 2. The current-voltage characteristics of the same annular JTJ as in Fig.l with 
one trapped fluxon and fluxon-antifluxon pairs. 



As a result, they leave a clear signature on the junction current-voltage 
characteristic (IVC). Fluxons having different topological charge n = ±1 
travel in opposite directions. Quantitatively, if a fluxon travels around an 
AJTJ having circumference C with a constant speed v < c(T), then an 
average voltage V = $>ov/C develops across the junction. By changing the 
bias current through the barrier the voltage drop changes and a new branch 
appears on the junction IVC. When N fluxons travel around an AJTJ, the 
junction voltage is V = N&qC/v. In this expression N is the total number 
of travelling fluxons and can be larger than the winding number n if FF 
pairs are travelling around the annulus. Therefore, we can count the number 
of travelling fluxons by simply measuring the voltage across the AJTJ. 

To show this, Figs.l and 2 represent the IVC of the same AJTJ with 
no fluxon trapped and with one fluxon trapped, respectively. We note 
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that with no trapped fluxons the zero voltage current is very large. In the 
other case the supercurrent is rather small (ideally zero) and large current 
branches can be observed at finite voltages corresponding to the fluxon and, 
possibly, FF pairs travelling around the junction. 

2.2. THE EXPERIMENTAL SETUP 

High quality Nb/Al — Al ox /Nb JTJs were fabricated on 0.5mm thick silicon 
substrates (chips) using the trilayer technique described in Ref.[22]. Each 
JTJ had a mean circumference C = 500 fj,m and a width Ar = 4 \xm. The 
thicknesses of the base, top and wiring layer were 200, 80 and 400 nm, 
respectively. The high quality of the samples was inferred from the I-V 
characteristic at T = 4.2 K by checking that the subgap current I sg at 2 mV 
was small compared to the current rise AI g in the quasiparticle current at 
the gap voltage V g . 

The symmetry of the junctions is assured by the absence of a logarithmic 
singularity in the IVCs at low voltages and the linear temperature depen- 
dence of the critical current as the temperature T approached the critical 
temperature Tq. Many samples have been measured. For clarity only two 
will be discussed here. Their geometrical and electrical (at 4.2 K) parame- 
ters of are listed in Table I. In particular, they differ in their critical current 
densities. The Swihart velocity cq has the value cq = 1.4 x 10 7 m/sec. 



TABLE 1. The samples 



Sample 


A 


B 


Mean circumference C(fj,m) 


500 


500 


Width Ar(^m) 


4 


4 


Zero field critical current I (mA) 


33 


2.5 


Maximum critical current / maJ< (myl) 


39 


2.7 


Gap quasiparticle current step AI g (mA) 


88 


5.2 


-^max / AIg 


0.45 


0.52 


Critical current density J c (A/cm 2 ) 


3050 


180 


Josephson length Xj(fim) 


6.9 


28 


Normalized mean circumference C/Xj 


72 


18 


Quality factor V m (mV) 


49 


63 


Normal resistance _Rjv(mf2) 


36 


610 


ZFS1 asymptotic voltage (fiV) 


51 


53 



In order to vary the quenching time in the broadest possible range, 
we have realized the experimental setup schematically shown in Fig.3. A 
massive Cu block held to the sample holder was used to increase the system 
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Figure 3. Sketch (not to scale) of the cryogenic insert developed to perform the junction 
thermal cycles over slow and fast timescales 



thermal capacity. The chip was mounted on one side of this block. On the 
other side, a thermoblock was mounted consisting of a 50 Q. carbon resistor 
and two thermometers in order to measure and to, if necessary, stabilize the 
Cu block temperature. Finally a small sized 100 S7 surface mounted resistor 
(SMR), was kept in good thermal contact with the chip. 

This system allowed us to perform the sample quenching over two quite 
different time scales. By means of the resistor in the thermoblock, a long 
time scale was achieved by heating the chip through the Cu block; on the 
contrary, a short current pulse through the SMR on the chip, attained 
much short thermal cycles. The total system was kept in a vacuum-tight 
can immersed in the LHe bath at He gas pipeline pressure. By changing 
the exchange gas pressure and using the two techniques, the quenching time 
can be changed over a quite large range from tenths to tens of seconds. 

The temperature dependence of the junction gap voltage was exploited 
to monitor the temperature of the junction itself during the thermal cycle. 
The analytical expression [20] for the superconductor gap energy A(T), 

A(T) , A(T)T C 

-A-f = tanh-AA-^, 4) 
A(0) A(0) T w 

also applies to the junction gap voltage that is proportional to it. In this 
way we derive temperature profiles for slow and fast cycles, assuming for 
A(0) and T c the values 2.85 meV and 8.95 K, respectively, as found by 
Monaco et al. [21] on similar JTJs. 
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Figure 4- Time dependence of the junction temperature (only reliable below the hori- 
zontal dashed line) during a slow thermal cycle, as determined from (4). The thick dashed 
line is the best fitting curve to (5) . 

We are only interested to the cooling process, and we use (4) to fit our 
data by a simple thermal relaxation equation: 

T(t) = T fin + (T tn -T fm )exp^- t -^^j (5) 

with only two fitting parameters to and r, and Tj n and Tfi n fixed at 8.95 
and 4.15 K, respectively. In Eq.5 to is the time at which T = Ti n = T c and 
r is the relaxation time which sets the cooling time scale. The fitting curve 
for a fast quench is shown by a thick dashed line in Fig.4, corresponding to 
a thermal relaxation time r of 0.073 s. The dashed horizontal line indicates 
the temperature threshold below which the temperature time dependence 
can be reliably accounted for by our measured data. The quenching time 
tq can be obtained from its definition tq = —Tq/T at T = Tc, giving 

TQ = T Tel (T in ~ Tf in ). 

2.3. THE MEASUREMENTS 

Quenching experiments were carried out in a double /U-metal shielded cryo- 
stat and the transitions from the normal to the superconducting states 
were performed with no current flowing in the heaters and the thermome- 
ters. Both the junction voltage and current leads were shorted during the 
all thermal cycle. Furthermore, the heat supplied to the sample was such 
that the maximum temperature reached by the junction was made slightly 
larger than its critical temperature, say at about 10 K, in order to make 
sure that also the bulk electrode critical temperature (Tc ~ 9.2K) was 
exceeded. In this case, according to Eq.5, the value of the quenching time 
is tq ~ 1.7 r. The value of the quenching times has been determined to an 



8 



0.2 



P 



0,1- 




0.05 



0,1 



1 



10 



x Q (sec) 



Figure 5. Log- log plot of the probability Pi to trap one fiuxon versus tq. The solid line 
is the best fitting curve, assuming a power law dependence as in (1). 

overall accuracy of better than 5%. For each value of the quenching time, 
in order to estimate the trapping probability, we have carried out a set 
of 300 thermal cycles and at the end of each cycle the junction IVC was 
inspected in order to ascertain the possible spontaneous trapping of one or 
more fluxons. 

As we shall see later, the AJTJs are such that the ZK causal length 
£ > C by roughly an order of magnitude. Thus the probability of finding a 
single fiuxon is small. In the following we will focus our attention only on the 
probability Pi to trap just one fiuxon, although a few times we found clear 
evidence of two and, more seldom, three homopolar fluxons spontaneously 
trapped during the N-S transition. However, these events were too rare to 
be statistically significant. Experimentally, we define P\ as the fraction of 
times in which at the end of the thermal cycle the junction IVC looks like 
Fig.lb, i.e. with a tiny critical current and a large first ZFS. Our definition 
of Pi is reasonable as far as the chance to trap two fluxons is negligible. 

3. The results 

When £ > C we estimate the probability of finding a fiuxon in a single 
quench to be 



where, from (3), a = 0.25. 

Fig. 5 shows on a log-log plot the measured probability Pi of a single 
fiuxon trapping obtained by quenching the sample A 300 times for each 
value of the quenching time tq changed by varying the exchange gas pres- 
sure and by using both the fast and slow quenching techniques. Although 




(6) 
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the points are quite scattered, we attempted to fit the data with an allo- 
metric function P\ = a ■ Tq b with a and b being free fitting parameters. 
We found that the best fitting curve, shown by the solid line in Fig. 6, has 
a slope b = 0.27 ± 0.05. Such a value of b, although affected by a 20% 
uncertainty, is in good agreement with the expected b = 0.25 dependence. 
For the coefficient a we found the best fitting value of 0.1 ± 10% (tq 

in seconds). This is to be compared with the predicted value of Ctq^/^o . 
Sample A had a circumference C = 500 /im. Its effective superconduc- 
tor thickness was d s 250 nm. At the final temperature Tf in = 4.2 K, 
the critical current density was J c (Tfi n ) = 3050 A /cm 2 and the Josephson 
length was Aj(Tjj ra ) = 6.9 fim. From this, and c$ given earlier, we infer 

that £o ~ 3.8 nm and tq rs 0A7ps. This then gives Ctq^/^q w 0.08 s 1 / 4 , 
in good agreement with the experimental value of a, given the fact that we 
only expect agreement in overall normalization to somewhat better than 
an order of magnitude level. 

Similar measurements have been carried out for sample B. Although 
compatible with (3), the results were affected by a data scattering even 
larger than that found for sample A. This is due to a much smaller normal- 
ized length which, according to (6), translates in a expected probability P\, 
for a given tq, about four times smaller, far too small to get statistically 
significant data in reasonable times. However the roughly measured prob- 
ability Pi of 1 fluxon every 50-100 attempts is in fairly good aggreement 
with the expected value. 

4. Comments, Future Experiments and Conclusions 

We consider this experiment to give a striking confirmation of the Zurek- 
Kibble predictions. 

Of the several experiments [7-13] that had been performed prior to ours 
to test (1) only two [12, 13] had variable tq and could confirm scaling 
behaviour. Of these, only one experiment [13] shows KZ scaling. Given this 
relatively poor success rate in confirming (3) we are considering a further 
experiment to measure the ZK scaling exponent, this time with manifestly 
non-symmetric AJTJs, which require different fabrication techniques. In 
this case the value of a inferred from the same causal arguments is a = 1/7, 
very different from the a = 1/4 that we tested above. 

Our experiments have demonstrated that quenching times of the order 
of 1 second give a rather large probability to trap one fluxons on AJTJs 
having a very large normalized length. However, very long junctions mean 
very large critical current densities that, in turn, require Josephson barriers 
so thin that their quality and uniformity is often spoiled. For these reasons, 
it would be highly desirable to compensate the reduced junction length 
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with an increased quenching rate, as it is suggested by the findings for 
sample B. On the contrary, in order to lower the quenching time, it is 
necessary to resort to new techniques since the maximum power that can 
be dissipated by the surface mount resistors sets an obvious lower threshold 
on tq. One possible way to reach this goal is to perform the junction thermal 

cycle by means of light pulses. Light dissipates inside the superconducting 
electrodes, but not in the substrate providing a local junction heating that 
will relax much faster to the background temperature. We estimate that, 
by using a properly focussed pulsed light beam, the quenching time scale 
can be reduced to the (is range. 
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